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ABSTRACT 

We studied several representative circumstellar disks surrounding the Herbig Ae star 
MWC 480 and the T Tauri stars LkCa 15 and DM Tau at (sub-)minimeter wavelengths 
m hues ofCClfl. 

Our aim is to characterize photochemistry in the heavily UV-irradiated 



^Based on observations carried out with the IRAM Plateau de Bure Interferometer. IRAM is supported by 
INSU/CNRS (France), MPG (Germany) and IGN (Spain). 
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MWC 480 disk and compare the results to the disks around cooler T Tauri stars. We 
detected and mapped CCH in these disks with the IRAM Plateau de Bure Interferome- 
ter in the C- and D-configurations in the (1-0) and (2-1) transitions. Using an iterative 
minimization technique, the CCH column densities and excitation conditions are con- 
strained. Very low excitation temperatures are derived for the T Tauri stars. These 
values are compared with the results of advanced chemical modeling, which is based 
on a steady-state flared disk structure with a vertical temperature gradient, and a gas- 
grain chemical network with surface reactions. Both model and observations suggest 
that CCH is a sensitive tracer of the X-ray and UV irradiation. The predicted radial 
dependency and source to source variations of CCH column densities qualitatively agree 
with the observed values, but the predicted column densities are too low by a factor of 
several. The chemical model fails to reproduce high concentrations of CCH in very cold 
disk midplane as derived from the observed low excitation condition for both the (1-0) 
and (2-1) transitions. 

Subject headings: stars: formation — planetary systems: protoplanetary disks - astro- 
chemistry - molecular processes - stars: individual (DM Tau, LkCa 15, MWC 480) 



Introduction 



Our present understanding of the chemical composition and evolution of protoplanetary disks 
is far from being complete. Recent years have revealed increased interest in such studies as a 
necessary prerequisite to characterize the physical and chemical conditions for planet formation. 
Apart from CO and its isotopologues, and occasionally HCO+, CN, HCN, and CS, the molecular 
content of prot oplanetary disks characterized by rnillime ter line observations remains poorly known 
(see reviews bv lDutrev et al.ll2007al : iBergin et al.l 120071 ). Recently, infrared spectroscopy provided 
evidence for the presence of C2H2, HCN, H2O, OH, and C O2 in the inner regions of protoplanetary 
disks and iridicated the importan . ce of photochemistry (e.g.. lLahuis et al.ll2006l : ICarr &: Naiitall2008l : 
Salvk et al.ll2008l : iPascucci et al.ll2009l ). Molecular line data are limited in their sensitivity and spa- 
tial resolution, which im plies that the spatial distribution of molecular abundances in disks remains 



poor ly determined (e.g., iPietu et al. 



2005 



Semenov et al.ll200,4 lOi et al.ll200fil . I2OO8I : IPanic et al 



2OO9I ). Thus a detailed comparison with existing chemical models is difficult and often based on 
global data (e.g. integrated line proffies), while for testing the validity of sophisticated models one 
has to include the thermal disk structure and the time-dependent chemical evolution together with 
line data and channel maps. 

As part of the Heidelberg-Bordeaux "Chemistry In Disks" (CID) project we are investigating 
the molecular content and physical properties of a sample of well-studied T Tauri and Herbig Ae 
disks of various age, followed by comprehensive physico-chemical modeling. In the first CID paper, 
we have presented the results of a deep search for N2H"^ and HCO"*" toward two T Tauri stars 
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(DM Tau, LkCa 15) and one Herbig Ae star (MWC 480), see lOutrev et alJ (|2007bl . CIDl paper 
hereafter). The N2H"'" emission has been detected in LkCa 15 and DM Tau, with the N2H+ to 
HCO"*" ratio of a few percent, similar to that of cold dense cores, and the disk ionization degree of 
~ 10~^, as predicted by chemical models. 



In the second CID paper by ISchreyer et al.l (j2008l ). we have found that the modeled and 
observed column densities of several key species relative to ^^CO are lower for the Herbig AB Aur 
disk than the values measured in DM Tau, while the absolute amount of CO gas is similar in both 
disks. This has been interpreted as an indication of a poor molecular content of the Herbig AOe 
system compared to the disk around the Mle star DM Tau due to intense UV irradiation of the 
disk by the AO star. 

In this third CID paper, we report on the observations of CCH in three objects: DM Tau, 
LkCa 15 and MWC 480. The CCH molecule has been selected because it should be sensitive to 
the UV radiation field and because its chemistry is relative ly well-studied, with many accurately 



determined relevant reaction ra tes (e.g., Woodall et al. 20071 ) . Using the ^^-minimization technique 



in the plane as described in iGuilloteau Sz DutrevI (j 19981 ). we derive the column density of CCH 
in the upper layers of the outer disks. We compare these values with the CCH column densities 
computed with realistic two-dimensional steady-state disk models and a gas-grain chemistry with 
surface reactions. 



2. Observations 

2.1. Sample of Stars 

The source properties and their coordinates are summarized in Table [TJ All these systems are 
isolated and located in regions devoid of CO contamination. The disk sizes, masses, and ages are 
presented in Table [2] (see also Tables 1 and 2 in CIDl). 



2.2. PdBI observations and data reduction 



The Plateau de Bure observations of the C2H(l-0) line in DM Tau were carried out in Septem- 
ber & October 1997, and between October 2002 and March 2003, mostly in compact configurations. 
The total observing time on source (equivalent to 6 antennas) is about 30 hours and the angular 
resolution (with natural weighting) is 3.5". The (1-0) line in LkCa 15 and MWC 480 was observed 
in track-sharing mode in December 2004 and April 2005, with a total observing time on source of 
~ 5 hours and an angular resolution of 3.8" x 2.5". The C2H(2-1) line was observed in track-sharing 
mode for the three sources in D configuration in August & September 2008. The total observing 
time on source with dual polarization receivers is 5 hours. The angular resolution is 2.8" x 2.5" 
(with natural weighting). 



Table 1. Properties of the targets 



Source 


Right Ascension 
(°' ,"), J2000.0 


Declination 
™,^), J2000.0 


Spect.Type 


Effective Temp. 
(K) 


Stellar Lum. 

(L©) 


Stellar Mass 
(M©) 


Age 
(Myr) 


UV flux 
(Xo) 


X-ray lum. 
ergs"-^ 


DM Tau 


04:33:48.73 


18:10:09.89 


Ml 


3720 


0.25 


0.55 ±0.07 


5 


300 


31029 


LkCa 15 


04:39:17.78 


22:21:03.34 


K5 


4350 


0.74 


0.97 ±0.03 


3-5 


1850 


31029 


MWC 480 


04:58:46.26 


29:50:36.87 


A4 


8460 


11.5 


1.65 ±0.07 


7 


5000 


10^8 



Note. — Coordinates J 2000.0 deduced from the fit of the 1.3mm PdBI cont inuum map llPietu et al.ll2006l ). Col. 3, 4, 5, 6 and 7 are the spectral type, 
effective temperature, t he stellar lumi nosity, the stellar mass, and age as given in lSimon et al. I ll2000h. T he ste llar FUV luminosities at r = 100 AU are given 
in Col. 8 in units of the lDrainj l ll978h interstellar UV field. These values are rescaled from lBerein et al.l l l2003h (Lk Ca 15 and DM Tau) for the 1 100 - 2 070A 
region, assuming that 30% of the FUV fiux is the Ly^ radiation. For MWC 480 the UV flux x is computed from the lKuriic3 l|l993l 'l ATLAS9 of stellar spectra. 
The stellar X-ray luminosities are taken from lCIassgold et al. I 1I2OO5I) . 
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Data were reduced with the GILDAS package. Flux caHbration was performed using MWC 349. 
The continuum was subtracted for the analysis. The channel width is 0.25 kms~^ for the J=l-0 
line and 0.26 kms~^ for the J=2-l line. Integrated spectra over a 7 x 7" region (for DM Tau) and 
5 X 5" region (for MWC 480 and LkCa 15) are presented in Fig. [TJ We detected C2H 1-0 emission 
from DM Tau and LkCal5, but not from MWC 480. The weak J=2-l transition is detected in the 
T Tauri sources (see Fig. [T]). The apparent disagreement between the observed and modeled line 
profiles for the DM Tau and LkCa 15 spectra are due to the relatively low signal-to-noise (SNR) 
level of our PdBI data. 

Channel maps of the C2H (1-0) emission in DM Tau are presented in Fig. [2l Since the signal- 
to- noise ratio of the C2H (1-0) spectra from MWC 48 and LkCa 15 is low , we analyzed these data 
as in our study of N2H~*" J=l-0 in the same sources (IDutrey et al.l l2007bl ) . A filter to the spectra 
is applied, in which each velocity channel is weighted by the expected line intensity (obtained from 
the best fit modeling), and the sum over the channels is normalized by the sum of the weights 
(a classical optimal filtering). The resulting SNR maps are presented in Fig. [3l which illustrate 



marginal detection of C2H (1-0) in LkCa 15 and MWC 480 (IDutrey et al 



2007 



2.3. Derived column densities 

The distribution of C2H in DM Tau, LkCa 15 and MWC 480 was analyzed using a parametric 
disk mo del, in which the column density and excitation temperatures are power laws of the ra- 



dius (see iPietu et al.l l2007l . for a more detailed description of the model). We used the geometric 
parameters (position, inclination, Keplerian velocity law ) derived from the CO isoto pologues and 
HCO^ (see Table [2]). The x^-minimization technique of Guilloteau &: Dutrey (jl998l ) was applied 



to constrain the CCH column densities and other disk parameters. The observation of two C2H 
lines allows us to recover the excitation conditions (within the limitations due to low signal to noise 
data). As the lines are relatively weak, the exponent of the excitation temperature is not well 
constrained. We find q = 0.4 ± 0.3 for DM Tau, and assumed q = 0.4 for the other sources. Such 
a value is in agreement with results obtained from stronger lines, in particular HCO^ (Table [2]). 
Consequently, the excitation temperature of the DM Tau disk is well constrained, and poorly de- 
rived for LkCa 15. For MWC 480 we had to assume T = 14 K in the data analysis. The obtained 
column densities are insensitive to the choice of the temperature in a range between 10 and ~ 20 K 
at 100 AU, and scale linearly with T for higher temperatures. The slope of the surface density 
distribution is constrained for DM Tau and LkCa 15. For MWC 480 we used a similar value. The 
derived temperature and column densities are summarized in Tabled 
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Table 2. Properties of protoplanetary disks 



Source 


i 


Rout 


Tioo 


q 


Acc. Rate 


Disk Mass 




n 


(AU) 


(K) 




(Moyr-i) 


(Mo) 


DM Tau 


-32 


800 


15 


0.12 - 0.5 


2.10-9 


0.05 


LkCa 15 


52 


680 


15 


0.4 


10-8 


0.03 


MWC 480 


38 


695 


30 


0.28 - 0.65 


10-8 


0.03 



Note. — Col. 2 Inclinations are taken from lSimon et al. (200C). Col. 3 & 4 Apparent di sk radii and 



tempe rature distribution are taken from the CO and HCO"*" interferometric analysis of lPietu et al 
(j2007l ). Here Tioo denotes gas kinetic temperature at 100 AU, q - the power-law exponent in the 



tempe rature distribution. Col. 5 Accretion rates are taken as a generic value from lD'Alessio et al 



1999 ) for MWC 480 and LkCa 1 5. The accret i on ra te for DM T au is taken from 



Calvet et al 



(120051). Col.6 Disk masses arefrom iDutrev et al.l (|l997l ) for DM Tau, iThi et al.l (120041 ) for LkCa 15 



and [Mannings k SargentI (| 19971 ) for MWC 480. 



Table 3. Derived column densities of CCH 



Source 


DM Tau 


LkCa 15 


MWC 480 


S at 300 AU (10^^ cm-^) 


2.8 ±0.2 


2.9 ± 1.1 


< 0.8 ±0.3 


P 


0.5 ±0.2 


0.5 ±0.2 


[0.5] 


Rout (AU) 


780 ± 40 


570 ± 50 


[400] 


T at 300 AU (K) 


7.3 ±0.6 


6.3 ±1.4 


[9] 


q 


0.4 ±0.3 


[0.4] 


[0.4] 



Note. — Values between square bracke ts indicate fixed p arameters as constrained by the previous 
CO, HCO+ and new CCH data fitting. dPietu et al.ll2007l ). The derived CCH column density for 
MWC 480 is an upper hmit. 
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Fig. 1.— Integrated spectrum of C2H J = 1 - (left) and J = 2 - 1 (right) for DM Tau, LkCa 15 
and MWC 480 integrated over the 7" x 7" central region. The thick curves represent the best-fit 
model deduced from the minimization procedure. 
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Fig. 2. — Channel maps of the C2H J = 1 — emission towards DM Tan. Contour spacing is 
11 mJy/beam, or 0.17 K {2a). Velocities are given in the upper left corners of the panels. The 
cross indicates the orientation of the disk. The beam sizes are indicated in the lower left corners. 




Fig. 3.— (Top to bottom) SNR maps of the C2H J = 1 - (left) and J = 2 
DM Tau, LkCa 15 and MWC 480. The contours correspond to the 2 a signal, 
indicated in the lower left corners. 



— 1 (right) toward 
. The beam size is 
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3. Disk physical and chemical model 



We refrain from comparing deconvolved observed and simulated spectral maps and base our 
analysis on the observed q uantities derived with the x^-™iiii™ization method in the uv-plane 
(jGuilloteau Dutrevlll998l ). This allows us to exclude computationally expensive iterative fitting 
of the data using 2D line radiative transfer modeling. 

For all sourc es, we adopted a flare d 1+lD-disk model with a vertical temperature gradient 
similar to that of lD'Alessio et alj (|l999l ). see Fig. [H The a-parameter for all three disks is chosen 
such that the modeled disk masses are in accordance with the observed values. Input disk and stellar 



parameters are either determined by the y -fitting or taken frorn prev i ous works ( Dutrey et al.lll997 : 



Mannin gs fc Sargent 



2005; Pietu et al. 



1997 



D'Alessio et al 



1999 



Simon et al. 



20001 : 



Thi et al. 



2004 



Calvet et al. 



20071). see Tabled The adopted viscosity parameter a is 0.01 for DM Tau, 0.03 



for LkCa 15, and 0.02 for MWC 480. It is chosen such that the resulting disk masses are the same 
as derived from millimeter dust continuum measurements. The dust grains are modeled as compact 
spheres of uniform 0.12 jum radius made of amorphous silicates with olivine composition and the 
optical data taken from lDraine fc: Lee ( 1984 ). The standard 1% dust-to-gas mass ratio is assumed. 
In addition, we consider a model of the DM Tau disk with bigger 1.0/im grains. 

The disks are illuminated by the UV radiation from the central star and by the inte rstellar 
UV ra diation. The intensity of the interstellar UV radiation field is x = 1 as derived by iDraine 
(j 19781 ) . The UV intensity at a given disk location is calculated as a sum of the stellar and interstellar 
components that are scaled down by the visual extinction in vertical direction and in direction to the 
central star (ID plane-parallel approximation). The non-thermal FUV radiation field from DM Tau 
and LkCa 15 is represented by the scaled ISRF, while that of MWC 480 is repre sented by the scaled 



radiat ion field of the Herbig AO star (according to Ivan Dishoeckl (|l988l ) and Ivan Dishoeck et al 



(j2006l )). We model the attenuation of cosmic rays (CRP) by Eq. (3) from ISemenov et al.l (|2004l ). 
In the disk interior, ionization due to the decay of short-living radionu clides is taken into account. 



assuming an ionization rate of 6.5 -10 s ([Finocchi &: Gaill Il997l ). The inte nsity and radial 



penet ration of the stellar X-ray radiation are m odeled using obs e rvation al results of lGlassgold et al, 
(|2005l ) and the 2D Monte Carlo simulations of lciassgold et al.l (!l997i 



We did not find any relevant data regarding observed X-ray luminosities for our stars and 
used the following estimates: Lx = 10^^ ergs~^ for MWC 480, and 310^^ ergs~^ for DM Tau and 
for LkCa 15. The low value of the adopted X-ray luminosity of MWC 480 reflects the fact that 
Herbig Ae stars have weak surface magnetic fields due to their non-convective interiors and thus 
lack an aoj dynamo generation mechanism. Consequently, thei r coronal activity is l ower and X-ray 
luminosities are 10-1000 times lower than that of T Tauri stars (jCiidel &: Nazell2009l ). Indeed, a few 
measurements of the strength of the longitudinal magneti c field in MWC 480 resulted in marginal 
detection of i? = 87 ± 22 Gauss (IHubrig et al.ll2004l . l2006l ). For the T Tauri sources we adopted a 
median value as recently measu red with Chandra and XMM in the Orio n Nebula Cluster and the 
Cepheus B star-forming region (jPreibisch et al.ll2005l : ICetman et al.ll2009l ): log(Lx/-Z^boi) ~ —3.5 = 
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310 ergs (with uncertainty of an order of magnitude). To investigate the role of the stellar 
X-ray radiation for the CCH chemical evolution, we consider three additional models of the DM 
Tau disk with Lx = 0, 10^9 and lO^o ergs'^ 

We added a new class of the X-ray driven reactions to our chemical code leading to the 
formation of doubly-charged atoms (O"^"'", C"'""'", N"*""*", S~^~^, Fe"^"'", Si" !""*") due to the Auge r effect . 
Secondary electron impact ionization cross sections were taken from iMeiierink &: SpaansI (|2005l ) . 
Neutralization reactions of thes e doubly-charged ioni zed atoms by e- and charge transfer reactions 
with molecules are adopt ed from lStauber et al.l (120051). The pho toionization cross sections are taken 
from lVerner et al.l (|l993l ). as described in lMalonev et al.l (|l996l ). In addition, we duplicated a few 
tens reactions of CRP with molecules to get the same reaction set for the X-ray photons. The rate 
of C O was modified in or der to take self-shielding properly into account according to the Eq. 9 
from lMaloney et al.l (jl996l ). 



The gas-grain time-dependent chemical model adopted in this study is mostly the same as 
in the CID2 paper. First, the osu.2007 network of gas-phase reaction^ is used. Second, we uti- 
lize a standard ra te approach to the surface chemistry modeling but without H and H2 tunneling 



(Katzetal 



model of 



19991). The sur f ace re actions together with desorption energies were a dopted from the 



Garrod &: HerbstI (120061). w hich is mostly based on previous studies of iHasegawa et al 



(jl992l ) and iHasegawa &: HerbstI (Il993l ) . Finally, rates of several tens of pho todissociation and pho- 
toionization reactions are updated according to Ivan Dishoeck et al.l (j2006l ). using the UV spectral 
shape typical of a Herbig Ae star and ISM-like dust grain optical properties. Overall, the disk 
chemical network consists of more than 650 species made of 13 elements and 7500 reactions. 



Using this time-dependent model and the "low metals" initial abundances of lLee et al.l (|l998l ). 
distributions of the molecular abundances and column densities for the considered species are 
simulated over 5 Myr of evolution. The choice of initial abundances does not affect the resulting 
molecular abundances and column densities due to relatively hig h densities in and long evolutionary 
timescales of protoplanetary disks (see also lWillacy et al.ll2006l ). 



4. Results of chemical modeling 

In Figure |5l the radial and vertical distribution of the absolute and relative C2ll abundances 
(to the total amount of hydrogen) at 5 Myr are shown. The highest relative C2H abundance is 
only about 10"^ (in DM Tau and LkCa 15) and several times lower in MWC 480. In DM Tau 
and LkCa 15 the warm layer of ethynyl is located at ~ 1 pressure scale height, with a local peak 
towards r = 10 AU. Note that there is a second peak, better visible in relative concentration, 
which is located in the outer disk region and at higher disk heights (r > 500 AU, z/Hj. ~ 1.5), 
where the IS UV photons can penetrate. CCH is partly frozen out in all T Tau disk midplanes 



^See: https://www.physics.ohio-state.edu/r-^eric/research.html 



- 12 - 




R, AU 

Fig. 4. — Distributions of the dust temperature (left panel) and particle density (right panel, 
logarithmic scale) for the DM Tau, LkCa 15, and MWC 480 disk models (top to bottom). The gas 
and dust temperatures are assumed to be the same. The Y-axis is given in units of the pressure 
scale height (as derived for the disk atmosphere). 
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as its evaporation temperature is ~ 50 K. In the warmer disk of MWC 480 ethynyl reaches the 
highest relative abundance only in the outer disk region, at ~ 1, while a peak in absolute CCH 
concentration is reached in the midplane around 100 AU. What is the reason for this chemical 
diversity of CCH in all 3 disks? 

In essence, the chemical evolution of the ethynyl radical is determined by removal rates of 
carbon and oxygen into more complex molecules, e.g. CO, OH, H2O, as well as by photoprocesses 
and freeze-out/evaporation, see Fig. [71 Our chemical simulations begin with atomic initial abun- 
dances. Thus, light ionized hydrocarbons CH^ (n=2..5) are the first molecular species produced 
by radiative association with H2, following hydrogen addition reactions: C^ CH^ CHj|" — )• 
CHg". The protonated methane reacts with electrons, CO, C, OH, and more complex species at 
later stage and forms methane. The CH4 molecules undergo reactive collisions with C"*", produc- 
ing C2H^ and C2H3'. An alternative way to produce C2H^ is the dissociative recombination of 
CH^ into CH3 followed by reactions with C^. The C2Hj^ ion is also formed by hydrogen addition 
to C2H^. Finally, C2H2' and C2H^ dissociatively recombine into CH, C2H, and C2H2. Another 
important formation channel for ethynyl is by surface hydrogenation of precursor species, C2. 

The major removal pathway for C2H is either the direct neutral-neutral reaction with O that 
forms CO, or the same reaction but with heavier carbon chain ions that are formed from ethynyl by 
subsequent insertion of carbon. In upper disk regions photodissociation of CCH into C2 and H (with 
the rate constant k = 5 .1 10~^° s~^ arid Ay exponent 7 = 1.9) and photoionization {k = lO"^'^ s~^, 



7 = 2.0) are efficient (Ivan Dishoeckl Il988l ). At later times, in colder disk regions depletion and 
gas-phase reactions with more complex species may enter into this cycle. One can restore high 
concentrations of CCH obtained at early times, t — 100 years by destroying CO and enriching 
the gas with elemental carbon. In disk it happens in the inner region either by reactions of CO 
with ionized helium atoms produced by the stellar X-ray radiation, or directly upon irradiation of 
CO by UV photons. This second rise in the CCH concentration is typically reached only at very 
late times in our models, t > 1 Myr. 

The relative importance of X-ray and UV radiation is shown in Fig. [6l using the disk around 
DM Tau as an example. The prominent intermediate layer of CCH almost disappears in the model 
without X-ray irradiation, and when grains grow beyond ~ 1/im. The lower opacities in the model 
with "big" grains shifts the upper outer layer of CCH towards midplane, following the enhanced 
penetration of the dissociating interstellar UV radiation. Thus CCH is indeed a sensitive tracer 
of the UV radiation field. The lack of X-ray irradiation results in a less "sharp" inner shape 
of the upper CCH layer (compare left and right panels in Fig. [6]), and low CCH abundances in 
the intermediate layer. In the absence of X-ray radiation, ionized helium is only produced by 
interactions with cosmic ray particles, and thus replenishment of C"*" and O into the gas upon 
destruction of CO by He^ is lower. Still, there is an "island" of high CCH concentration in the 
intermediate layer at r ~ 200 — 300 AU in the model of DM Tau without X-rays. 



MWC 480, as a Herbig Ae star, emits ~ 3 and ~ 15 times more UV photons than LkCa 15 
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and DM Tau, respectively, but has a weak X-ray luminosity (see Tabled]). Consequently, the upper 
CCH layer is more noticeable in the MWC 480 disk, where a larger fraction of elemental carbon can 
be kept from being locked in CO and other species by photodissociation. Note that the molecular- 
poor atmosphere is slightly more extended in the MWC 480 disk compared to those of LkCa 15 
and DM Tau (upper panel in Fig. [5|). The low X-ray luminosity of MWC 480 assumed in our 
modeling also leads to lower concentration of CCH in the inner warm molecular layer (r < 50 AU, 
~ 1 pressure scale height). However, there is a region around the midplane (50 — 200 AU, up to 
z/Hr ~ 1) where absolute abundances of CCH are high in the MWC 480 disk. Since the Herbig 
Ae disk is warm, surface processes involving heavy radicals are more active there, while volatile 
ices are nearly absent. Using our chemical analyzing software, we found that the oasis of CCH 
in MWC 480 is due to efficient surface hydrogenation of C2 ice, and other surface reactions (e.g., 
photodissociation of C2H2 ice, and surface formation of C2 from two C atoms). Temperature there 
ranges from 23-50 K, while densities are high, nn ~ 10^ — 10^ cm^^, resulting in fast accretion 
rates of precursor molecules to grain surfaces and sufficiently high desorption rate of CCH. Note 
that this oasis of ethynyl appears only after ~ 1 Myr of evolution, due to slow surface chemistry. In 
both T Tau disks photodissociation of CO is effective in the outer uppe r disk region. Both DM Tau 



and LkCa 15 have non- negligible X-ray luminosity of ~ 10^^ ergs ^ ( Glassgold et al.ll2005l ). and 
thus CCH can be re-formed toward inner disk regions, see Fig. [5l 

The observed and modeled CCH column densities are compared in Fig. [8] (left panel). The 
CCH column densities in the DM Tau disk computed with the model without X-ray radiation and 
the model with bigger l^m. grains are depicted in Fig. [8] (middle panel). The influence of the 
adopted stellar X-ray luminosity on the CCH column densities in the DM Tau disk is shown in 
the right panel. As we mentioned above, the X-ray photons enhance CCH abundances, leading to 
higher CCH column densities (by a factor of 3-10). In the model with micron-sized grains CCH has 
high column density only in the very outer disk region, at radii larger than ~ 300 AU, while it drops 
quickly in the inner disk zone. Given the uncertainty of the modeled column densities (a factor 



of ~ 3, see lWakelam et al.ll2005l : IVasvunin et al.ll2008l ) and observational uncertainties (see error 



bars in Fig. [8|), our relatively simple steady-state disk model is in a moderate agreement with the 
observed quantities. In particular, the flat radial slope of the CCH column densities is reproduced 
by the simulations. Note that the column density of CCH in MWC 480 is lower compared to that 
of the DM Tau and LkCa 15 disks, and a similar trend exists in the observational data. 

The effect of the stellar X-ray radiation field on the computed CCH column densities is such 
that the lower the adopted X-ray luminosity, the lower the resulting iV(CCH), see Fig. [HI right 
panel. Apparently, a transition from Lx = 10'^'^ ergs^^ to 3 lO'^^ ergs~^ is less prominent in the 
CCH column densities than a transition from Lx = 3 lO'^^ ergs~^ to lO'^^ ergs~^. Thus, indeed the 
ethynyl radical can serve as a sensitive probe of the impinging X-ray radiation field. In general, it 
is the stellar UV luminosity that affects the evolution of CCH globally, with X-ray radiation being 
only important in inner disk regions (r < 50 — 200 AU). 
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Discussion 



As we have shown in lSchrever et al.l (|2008l ). the stehar UV luminosity can be a crucial factor 
that makes the AB Aur system chemically poorer with respect to the disk of DM Tau. The same 
trend is found here for C2H and the disks around MWC 480 and DM Tau. However, we also show 
that X-ray radiation plays an important role for high CCH concentration in the inner regions of 
T Tau disks. Furthermore, the results of the steady-state disk chemical model is in qualitative 
agreement with the CCH column densitie s obtained by x^-fitting of the spectral maps. In contrast 
to our CIDl paper (IDutrey et al.ll2007bl ). where such a model could not fit well the radial slope 
of the N2H"'" column densities, in the case of CCH the flat distribution is obtained both from the 
analysis of the interferometric data and chemical modeling. 

Our chemical modeling predicts that two CCH molecular layers are located at elevated disk 
heights, or in warm disk regions (T ~ 35 — 75 K). This is inconsistent with the very low excitation 
temperature of the CCH emission derived for DM Tau. The critical densities for the CCH (1-0) and 
CCH(2-1) lines to be thermalized are 210^ cm~^ and 710^ cm~^, respectively. The two layers of 
CCH have higher densities of ~ 10^ — 10* cm~^ in all studied sources (compare Fig. [Hand Fig. [5]). 
Note that CCH layers are at higher densities in the MWC 480 disk compared to those around 
LkCa 15 and DM Tau. Hence, according to our model CCH(l-O) and (2-1) should be thermalized 
in all considered sources, and sub-thermal excitation of these lines is unlikely. 

Which factors that are not taken into account in our simplified chemical modeling could lead 
to a substantial concentration o f CCH in the cold disk midpl ane of DM Ta u and LkCa 15? As 
we have found in lSemenov et al.l (120061 ) (and later confirmed by I Aikawal 120071 ) . turbulent transport 
can enhance the CO concentration in the DM Tau disk midpla. ne such that ^^CO (1-0) and (2-1) 
lines are indeed "cold" (Tex ^ 15 K). Later. iHersant et al.l (|2009l ) have shown that the scattering of 
the stellar UV radiation by PAHs and dust grains in upper disk atmosphere can redirect a portion 
of these UV photons into the disk, allowing CO to efficiently desorb from dust grains in the cold 
disk zone. Despite the fact that this assumption is not considered in our model, where desorption 
of surface molecules in the disk midplane is only due to the CRP-driven UV radiation and CRP 
heating, a non-negligible amount of CCH exists in the outer disk midplane at T < 20 K, see Fig. [9] 
(left panel). The CCH column density in DM Tau at T < 20 K has a maximum value of 6 10^^ cm~'^ 
at r = 800 AU, which rapidly decreases with radius. In warmer regions of the DM Tau disk the 
CCH column density is about 5 times higher, ss 3 10^^ cm~'^, and has an almost flat profile. Thus, 
we expect that only < 20% of the CCH (1-0) and (2-1) line fluxes arise from the cold disk zone. 



The amount of cold gases in disk midplanes can also be increased if gas-gr ain interaction s are 
ineffective such that freeze-out is delayed toward late evolutionary times (e.g.. I Aikawal 120071 ) . We 



did several test calculations with bigger grains of l//m and found that indeed neutral species become 
more abundant in cold disk region. The results for molecular chemistry in disks using a realistic 
grain evolutionary model will be reported in a separate publication. In essence, in evolving disks 
substantial amount of gas-phase molecules in very cold region can be produced and sustained by 
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moderate grain growth and sedimentation of micron-sized and bigger grains toward the midplane, 
while a fraction of smaher grains and PAHs remai ns in the disk atmospher e and provides enough 
opacities to keep t he disk interior coo l (see e.g., lAikawa &: Nomural l2006l : iJonkheid et al.l 120071 : 
Brauer et allbood : birnstiel et aLlbood l. 



Finally, in our disk model dust thermal balance is used to compute the temperature distribu- 
tion, which is assumed to be the same for the gas. In the disk atmosphere gas will be hotter due 
to e.g. photoelectric heating of small dust grains/PAHs by stellar/IS UV &: X-rays photons. We 
do not take this effect into account. For CCH lines it is of no concern as derived low excitation 
temperatures for the 1-0 and 2-1 transitions are indicative of a substantial amount of ethynyl being 
located in the cold disk zone around the midplane. In addition, our modeling results show that the 
majority of CCH is located at scale heights of < 1 for DM Tau and LkCa 15 and is even closer to 
the disk midplane for MWC 480 (Fig. \5\ bottom row). When compared to the corresponding disk 
density structures (Fig. it becomes apparent that the CCH molecular layers are within the disk 
regions having particle densities of several 10^-10^ cm . These regions are dense enough to allow 
effective coupling of gas and dust kinetic temperatures. 



6. Summary & Conclusions 

We observed protoplanetary disks around the Herbig Ae star MWC 480 and the T Tauri stars 
LkCa 15 and DM Tau in the (1-0) and (2-1) lines of CCH. We detected and mapped these disks 
in the CCH lines, using the IRAM Plateau de Bure Interferometer in the C- and D-configurations. 
Using the iterative minimization technique, the column densities of ethynyl are derived, which are 
lower in MWC 480 (A^(CCH) < lO^^ ^,^-2^ ^^ian in DM Tau and LkCa 15 (iV(CCH) ~ 3 10^^ cm^^). 
The derived excitation temperature of CCH in DM Tau is very low, Tex ~ 6 K, which is along with 
previous detections of cold gas in this object (CO, HCO'^, CH, and HCN). The observed CCH 
column densities are compared with the results of advanced chemical modeling, which is based on 
a steady-state flared disk structure with a vertical temperature gradient, and a gas- grain chemical 
network with surface reactions. The CCH abundances in all disk models show two layers of high 
concentration. Overall, the disk around the Herbig Ae star MWC 480 has less CCH, which we 
explain by the strong UV radiation field of this object, and the lack of X-ray irradiation from the 
star. The modeled CCH column densities are in qualitative agreement with the observed values at 
an evolutionary time of a few million years. Altogether, our measurements and modeling support 
the observational trend that disks around Herbig Ae stars tend to be more molecule deficient 
compared to those surrounding T Tauri stars due to stronger UV irradiation and lower X-ray 
luminosities (affecting inner disk chemistry). Our model reproduces the radial slope of the CCH 
column densities. Yet it fails to explain the low excitation temperature of the observed CCH lines 
as the modeled CCH molecular layers are located in disk regions dense enough to thermalize both 
transitions. The chemical simulations show that only < 20% of the CCH column densities are in the 
region of the DM Tau disk that is colder than 20 K. Other mechanisms (stronger photodesorption. 
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Fig. 5. — (Top row; left to right) The CCH abundances (relative to the total amount of hydrogen 
nuclei) in the DM Tau, LkCa 15, and MWC 480 disks at 5 Myr. The vertical axis is expressed in 
units of a pressure scale, H{r) = ^/{2)cs/^, where Cg is the sound speed in the disk atmosphere 
and is Keplerian angular velocity. (Bottom row) The same as in the top row but for the absolute 
CCH concentrations (cm~'^). 




Fig. 6. — (Left to right) The relative CCH abundances in the DM Tau disk computed for 3 distinct 
models: (1) the standard model (top left panel, Fig. ED, (2) the model with bigger l^m grains, and 
(3) the model without X-rays. 
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Fig. 7. — Most important chemical pathways for the evolution of CCH in protoplanetary disks. 
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Fig. 8. — (Left) Radial distribution of the CCH vertical column densities in the disks of DM 
Tau (solid line), LkCa 15 (dotted line), and MWC 480 (dashed line) at 5 Myr. The observed 
column densities and their uncertainties are indicated by the lines with error bars. (Middle) Radial 
distribution of the CCH vertical column densities in the DM Tau disk computed with: (1) the 
standard model (solid line), (2) the model with bigger 1/im grains (dotted line), and (3) the model 
without X-rays (dashed line). (Right) Radial distribution of the CCH vertical column densities in 
the DM Tau disk computed with: (1) the Lx = 310^^ ergs~^ model (solid line), (2) the Lx = 
10^^ ergs""^ model (dotted line), and (3) the Lx = lO^'^ ergs~^ model (dashed line). 
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dynamic processes, grain growth) are required to enhance the abundance of CCH in the cold disk 
midplanes. 
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Fig. 9. — Radial distribution of the CCH vertical column densities at 5 Myr in the disk regions 
with T < 20 K (left panel) and T > 20 K (right panel) for DM Tau (solid line), LkCa 15 (dotted 
Une), and MWC 480 (dashed line). 



